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Abstract 



We present a new method to account for coherence length effects in a semi- 
classical transport model. This allows us to describe photo- and electropro- 
duction at large nuclei (A > 12) and high energies using a realistic coupled 
channel description of the final state interactions that goes beyond simple 
Glauber theory. We show that the purely absorptive treatment of the fi- 
nal state interactions can lead to wrong estimates of color transparency and 
formation time effects in particle production. As an example, we discuss ex- 
clusive p° photoproduction on Pb at a photon energy of 7 GeV as well as K + 
production in the photon energy range 1-7 GeV. 

PACS numbers: 25.20.Lj, 24.10.Eq, 24.10.Nz, 25.75.Dw 
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I. INTRODUCTION 



In a high energy collision between two hadrons or a photon and a hadron it takes a finite 
amount of time for the reaction products to evolve to physical particles. During the collision 
process some momentum transfer between the hadrons or some hard scattering between two 
of the hadrons' constituents leads to the excitation of hadronic strings. The time that is 
needed for the creation and fission of these strings as well as for the hadronization of the 
string fragments cannot be calculated within perturbative QCD because the hadronization 
process involves small momentum transfers of typically only a few hundred MeV. One can 
perform an estimate of the formation time Tf in the rest frame of the hadron. It should 
be of the order of the time that the quark- ant iquark (quark-diquark) pair needs to reach a 
separation, that is of the size of the produced hadron (r^ 0.6 — 0.8 fm): 



During their evolution to physical hadrons the reaction products will react with reduced 
cross sections. This is motivated by means of color transparency: the strings and the sub- 
strings created during the fragmentation are in a color singlet state and therefore react with 
a cross section that increases with their transverse size. As a consequence the produced 
hadrons travel inside the nuclear medium with a reduced scattering probability during their 
formation time. Hence the formation time plays an important role in the dynamics of nuclear 
reactions, e.g. heavy ion collisions, proton and pion induced reactions as well as photon and 
electron induced reactions on nuclei. 

The latter two are of special interest because they are less complex than heavy ion 
collisions and, in contrast to hadron induced reactions, the primary reaction does in general 
not only take place at the surface of the nucleus but also at larger densities. Experiments at 
TJNAF and DESY, for example, deal with exclusive and semi-inclusive meson photo- and 
electroproduction at high energies. Large photon energies E 1 are of special interest because 
of time dilatation the formation length If in the rest frame of the nucleus can exceed nuclear 
dimensions: 



If one chooses the formation time to be Tf = 0.8 fm/c, the formation length in the rest frame 
of the nucleus will be about 30 fm for a 5 GeV pion and about 7 fm for a 5 GeV kaon or a 
7 GeV p meson. These lengths have to be compared with the typical size of nuclear radii, e.g. 
2.7 fm for 12 C and 7.1 fm for 208 Pb. Because it suppresses the final state interactions, the 
formation time has a big effect on photonuclear production cross sections at high energies. 
Turning this argument around, exclusive and semi-inclusive photoproduction of mesons on 
nuclei [1] offer a possibility to study these formation times if the FSI are well under control. 

Usually, the FSI are modeled within Glauber theory [2,3] and treated purely absorptive. 
A more realistic description of the FSI has to take also regeneration of the mesons studied 
into account. For this we use a coupled channel semi-classical transport model based on the 
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Boltzmann-Uehling-Uhlenbeck (BUU) equation. Originally developed to describe heavy ion 
collisions [4] at SIS energies it has been extended in later works to investigate also inclusive 
particle production in heavy ion collisions up to 200 AGeV and n [5] and p induced as well 
as photon and electron induced reactions in the resonance region [6]. Inclusive photoproduc- 
tion of mesons at energies between 1 and 7 GeV has been investigated in [7]. An attractive 
feature of this model is its capability to describe a large variety of very different reaction 
types in a consistent way. 

Since photon induced reactions are known to be shadowed (ct 7 a < Aa lN ) above photon 
energies of approximately 1 GeV [8-10], one needs a way to account for this shadowing 
effect in photoproduction. This is straightforward within Glauber theory, but it is not clear 
how to account for this initial state coherence length effect in a semi-classical transport 
model for the FSI. A first attempt at combining the quantum mechanical coherence in the 
entrance channel with our incoherent treatment of the FSI has been made in [7]. In Sec. II 
we present a new, improved way to implement shadowing in our model. In Sec. Ill the re- 
sults for exclusive p° photoproduction on 208 Pb at 7 GeV are compared with the predictions 
from simple Glauber theory. In Sec. Ill we also demonstrate the importance of the coupled 
channel treatment of the FSI by looking at the nuclear K + photoproduction cross section. 
We summarize our results in Sec. IV. 



In this section we describe how we account for shadowing in our model and only sketch 
the principles of the transport model itself. For a detailed description of the transport model 
used we refer to [11]. In our model the incoherent reaction of a high energy photon with a nu- 
cleus takes place in two steps. In the first step the photon reacts with one nucleon inside the 
nucleus (impulse approximation) and produces some final state X. In this process nuclear 
effects like Fermi motion, binding energies and Pauli blocking of the final state nucleons are 
taken into account. In the second step the final state X is propagated within the transport 
model. Except for the elastic vector meson and exclusive strangeness production (see [7]) 
we use the Lund string model FRITIOF [12] to describe high energy photoproduction on 
the nucleon. The same string model is also used to deal with high energy particle collisions 
in the FSI. Since FRITIOF does not accept photons as incoming particles, p dominance was 
used in [7] and the photon was passed as a massless p°. This led to an excellent description 
of charged particle multiplicities in 7p collisions. 

Here we generalize to vector meson dominance (VMD) [13] by writing the incoming 
photon state as 



and pass the photon as a massless p°, u or with a probability corresponding to the strength 
of the vector meson coupling to the photon times its nucleonic cross section cjvn 



II. MODEL 
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These probabilities follow from using (2.1) in the optical theorem: 
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As can be seen from (2.1) and (2.3) there is also a finite probability 



p( 7o ) = i- Yl p{ y) 



(2.4) 



v=p,u>,4> 



that the 'bare' photon 7 has to be passed to FRITIOF; the component 70 does not get 
shadowed in the nucleus. In a generalized VMD model for example the 70 includes contri- 
butions from heavy intermediate hadronic states which, as we show later, are not shadowed 
because of the large momentum transfer qy that is needed to put them on their mass shell. 
Since FRITIOF does not accept a 'bare' photon as input we replace it again by a vector 
meson V(—p°, u or 0), with the probability 



The particle production in FRITIOF can be decomposed into two parts. First there is 
a momentum transfer taking place between the two incoming hadrons leaving two excited 
strings with the quantum numbers of the initial hadrons. After that the two strings fragment 
into the observed particles. As a formation time we use 0.8 fm/c in the rest frame of each 
hadron; during this time the hadrons do not interact with the rest of the system. 

Up to now we did not take shadowing into account. In Glauber theory this corresponds 
to using only the left amplitude in Fig. 1, where a photon directly produces some hadron X 
at nucleon j. The left amplitude alone leads to the unshadowed incoherent photoproduction 
cross section 



where n{r) denotes the nucleon number density, h{f) the number density of nucleons with 
the correct charge to produce the hadron X and o X n is the total XN cross section. The 
exponential damping factor in (2.6) describes the absorption of the particle X on its way 
out of the nucleus. 

In Glauber theory shadowing arises from the interference of the left amplitude in Fig. 1 
with the second amplitude of order a em which is shown on the right hand side. In this process 
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the photon produces a vector meson V on nucleon i without excitation of the nucleus. This 
vector meson then scatters at fixed impact parameter b (eikonal approximation) through 
the nucleus to nucleon j and produces the final state meson X, leaving the nucleus in the 
same excited state as in the direct process. Off-diagonal scattering, where a vector meson 
V scatters into a different vector meson V', is usually neglected. In view of the upcoming 
discussion of exclusive p° photoproduction on Pb we state here the Glauber result [3] for 
the incoherent vector meson production cross section: 

/j2i f°° i I? \ ~ a VN J z °° dzn(b,z) 
do dzjn{b, Zjje z * 
J —CO 

dzin(b, Zi)— — (1 — iay)e tqyZi exp — -ay N (\ — iay) I dzn(b,z) 
-oo 2 \_ 2 J z . 

Again n(r) denotes the nucleon number density, <7vn the total VN cross section and ay the 
ratio of real to imaginary part of the VN forward scattering amplitude. The momentum 
transfer 

qv * if (2 ' 8) 

arises from putting the vector meson on its mass shell. A large momentum transfer qy is 
suppressed by the elastic nuclear formfactor and leads to less shadowing due to the oscillat- 
ing factor exp[iq v Zi\ in the integrand of Eq. (2.7). Note that q v is just the inverse of the 
coherence length, i.e. the distance, given by the uncertainty principle, which the photon 
can travel as a V fluctuation. For the quantities in (2.7) we use the parameterizations of 
Model I of Ref. [13] with which one obtains a very good description of the shadowing effect 
in nuclear photoabsorption down to the onset region [9,10]. 

One clearly sees from (2.7) how the FSI separate from the 'initial state interactions' of 
the photon. We now replace the exponential damping factor exp[— <Jvn dzn(b, z)], which 
corresponds to purely absorptive FSI, by a transport model. This allows us to incorporate 
a wider class of FSI. In addition we want to include events where the final vector meson V 
is not produced in the primary reaction but via sidefeeding. This leads to less reduction of 
the nuclear production cross section even with the same avN- We thus also need to know 
how all the other possible primary reactions are shadowed. We, therefore, start from (2.1) 
and use Glauber theory to calculate how the single V components of the photon change due 
to multiple scattering on the way to nucleon j where the state X is produced [3]: 

W j ))=(l- £ ^r)l7o>+ £ f (l-IVft))|V). (2.9) 

Here the (photon energy dependent) nuclear profile functions for the different vector meson 
components V are given as 



v 



dzifiib, Zi)-^-(l - ia v )e iqv ^~ Zj) exp --oyjv(l - iay) / dzn(b, z) 
■oo 2 2 J 
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Note that the 70 component is, by definition, not modified due to the presence of the nuclear 
medium. The cross section for the photon to react with nucleon j at position fj inside the 
nucleus can be deduced via (2.9) from the optical theorem: 



3)=(i- £ ^r) £ (f) 2 |i- r v(^)|V^. (2.11) 



Like for the photon in vacuum each term gives the relative weight for the corresponding 
photon component to be passed to FRITIOF. When integrated over the whole nucleus one 
gets from Eq. (2.11) the total incoherent photonuclear cross section 



dhjn{fj)a lN {fj) (2-12) 



which is shown in Fig. 2 by the dashed line together with the total nuclear photoabsorption 
cross section as derived from the optical theorem in [9]: 

a t °X = Aa lN - ^2 7TT I — I ~& e \ &vn( 1 - a v) 2 / d 2 b / dzi / dzjn(b, Zi)n(b, zj) 



xe iqv{zi ~ Zj) exp 



1 T J 

— cvat(1 — «o;y) / dz'n(b,z' 

2 iz t 



(2.13) 



Here fc and £v denote the momentum of the photon and the vector meson respectively and 
we again use the parameterizations of Model I of Ref. [13]. More than 90% of the difference 
between those two cross sections stems from coherent p° photoproduction. 

In Fig. 3 we show how strongly the p° and the <fi components of a real 20 GeV photon 
are separately shadowed in Pb. We plot the number density of the nucleons reacting with 
the V components of the photon 



ffifi) = n(fj)^- (-f) |1 - r v (f ] )\ 2 a VN (2.14) 

V-yN \9V J 



as a function of fj. One clearly sees that due to its smaller nucleonic cross section the 
component is less shadowed than the p° component at the backside of the nucleus. As a 
consequence strangeness production (e.g. K photoproduction), where the primary reaction 
is preferably triggered by the component of the photon, is less shadowed than, e. g., 7r 
photoproduction. This dependence of the strength of shadowing on the reaction type is 
new compared to the treatment of shadowing in [7] and can also be seen directly from the 
second amplitude in Fig. 1 because of the occurrence of the scattering process VN — > XN 
at nucleon j. 

As already mentioned above the purely absorptive FSI of the Glauber model are very 
different from the coupled channel description of a transport model. The transport model 
we use is based on the BUU equation that describes the time evolution of the phase space 
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density fi(r,p, t) of particles of type i that can interact via binary reactions [11]. In our case 
these particles are the nucleons of the target nucleus as well as the baryonic resonances and 
mesons (n, 77, p, K, ...) that can either be produced in the primary reaction or during 
the FSI. For particles of type i the BUU equation looks as follows: 

(0 dH OH \ 
Ft + -WW WW) mm = - fi - jM] - (2 ' 15) 

In the case of baryons the Hamilton function H includes a mean field potential which in our 
model depends on the particle position and momentum. The collision integral on the right 
hand side accounts for the creation and annihilation of particles of type i in a collision as 
well as elastic scattering from one position in phase space into another. For fermions Pauli 
blocking is taken into account in I coll via blocking factors. For each particle type i such 
a BUU equation exists; all are coupled via the mean field and the collision integral. This 
leads to a system of coupled differential-integral equations which we solve via a test particle 
ansatz for the phase space density (for details see [11]). Since the collision integral also 
accounts for particle creation in a collision the observed outgoing particle X cannot only be 
produced in the primary reaction, but can also be created by sidefeeding in which a particle 
Y is created first which propagates and then, by FSI, produces X. In addition the state X 
might get absorbed on its way out of the nucleus but be fed in again in a later interaction. 
Both cases can a priori not be ignored, but are usually neglected in Glauber models. 



III. RESULTS 

Exclusive vector meson photo- and electroproduction on nuclei is an ideal tool to study 
effects of the coherence length, formation time and color transparency. Exclusive p° elec- 
troproduction has been investigated in the HERMES experiment [14] at photon energies 
between 10 GeV and 20 GeV and virtuality Q 2 < 5 GeV 2 . The calculations for meson 
production on nuclei are usually done within Glauber theory [15]. As already mentioned 
above the FSI in Glauber theory are usually purely absorptive. This means that for the 
reaction 7A — > p°A* the primary reaction has to be ^yN — * p°N. If one treats the FSI via 
an absorptive optical potential one gets an exponential damping ~ exp[— a p ^ dzn(b, z)\ 
of the nuclear production cross section. 

We are presently working at incorporating also incoming virtual photons into the for- 
malism developed in Sec. II and at enlarging the configuration space for the FSI. Here we 
discuss, therefore, only results obtained with real photons and at a lower energy. In Fig. 4 
we show the results of our model for the mass differential cross section of incoherent, ex- 
clusive p° photoproduction on 208 Pb for E 1 = 7 GeV. In this case 'exclusive' means that 
the final state consists of a 7r + 7r~ pair and 208 bound nucleons. The solid line represents 
a calculation with the primary reaction restricted to 'yN — > p°N. It already includes the 
effects of shadowing, Fermi motion, Pauli blocking and the nucleon potential, but no FSI. 
The dotted line shows the effect of the FSI without a formation time of the p° in <yN -> p°N. 
The Glauber model yields quantitatively the same effect of the FSI. This means that in this 
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case (only 7iV — > p°N as primary reaction) FSI processes like p°N — + 7riV, 7riV — >■ p°iV, 
where the primary p° gets absorbed first and is fed into the outgoing channel by a later FSI, 
are negligible. If one assumes a formation time of tj = 0.8 fm/c for the p°, one gets the 
result indicated by the dash-dotted line. Due to the finite formation time there is consid- 
erably less absorption and the nuclear production cross section increases. If the observed 
spectrum looked like this, one would in Glauber theory be lead to the conclusion of a finite 
p° formation time . 

However, one will get a similar result with Tj = if one allows for other primary reaction 
besides 7^ — > p°N and uses a coupled channel model. This can be seen by looking at the 
dashed line in Fig. 4. We find that about 60% of the additional p° stem from inelastic p° 
production in the primary reaction, e.g. ^yN — > p°nN where the 7r gets absorbed during the 
FSI. We now apply an exclusivity measure like the one used in the HERMES experiment [14] 

2 2 

-2 GeV < AE = Py ~ mN < 0.6 GeV, (3.1) 
2m N ' v ; 

where is the nucleon mass and 

Py=Pn+Pj-P p (3.2) 

is the 4-momentum of the undetected final state. Here p 7 and p p denote the 4-momenta of 
the incoming photon and the detected 7r + 7r~ pair and Pn is the 4-momentum of the struck 
nucleon which, for the calculation of py, is assumed to be at rest. Using (3.1) leads to a 
decrease of the cross section (dash-dot-dotted line) because some of the inelastic primary 
events are excluded. If the exclusivity measure was good enough to single out only the elas- 
tic 7iV — * p°N reactions from the primary events, the dash-dot-dotted curve would coincide 
with the dotted line and Glauber theory would be applicable. Since this is not the case, one 
still extracts a too large formation time when using Glauber theory. 

One therefore needs a further constraint in addition to the exclusivity measure (3.1) 
which becomes apparent by looking at the differential cross section ^ in Fig. 5; the mean- 
ing of the lines is as before. For \t\ > 0.1 GeV 2 the full calculation with exclusivity measure 
(dash-dot-dotted line) gives the same result as the one with the primary reaction — > p°N 
and FSI (dotted line). In this kinematic regime Glauber theory can therefore be used. For 
|i| < 0.1 GeV 2 , however, the exclusivity measure (3.1) cannot distinguish between elastic p° 
photoproduction (jN — > p°N) and other primary reactions, e.g. inelastic p° photoproduc- 
tion (jN — > p°X, I ^ JV). At low values of |i| there exist many states X with invariant 
masses which are not excluded by the exclusivity measure (3.1). In addition we find that 
about 25% of the finally accepted p° in this t region are not produced in the primary reac- 
tion but stem from sidefeeding in the FSI. In the HERMES experiment one makes a lower 
\t\ cut to get rid of the coherent p° photoproduction contribution. In the case of lead and 
E 1 = 7 GeV the coherent part can be neglected above \t\ ~ 0.05 GeV 2 . For Glauber theory 
to be reliable one has to increase this threshold to approximately \t\ = 0.1 GeV 2 . Glauber 
theory can thus be trusted only under certain kinematic constraints. 
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So far we have discussed a case in which a strongly interacting particle, the p° meson, is 
produced. In the following we will now discuss the special effects that appear in a coupled 
channel treatment of the FSI when a weakly interacting particle, such as the K + meson, is 
considered. In Fig. 6 we therefore show the cross section for the reaction — > K + X in the 
photon energy range 1-7 GeV for 208 Pb. This reaction had already been investigated in [7]; 
the new treatment of shadowing and the initial interactions as outlined in Eqs. (2.1)-(2.5) 
lead to an increase of the total yield at 7 GeV by about 20%. The solid curve in Fig. 6 
represents the results of the full calculation (including shadowing, FSI, 77 = 0.8 fm/c, etc.). 
By comparison with the calculation without shadowing (dash-dotted line) one sees how im- 
portant shadowing becomes with increasing photon energy. At 7 GeV it reduces the nuclear 
production cross section to about 65%. 

The importance of a coupled channel treatment of the FSI becomes clear when compar- 
ing the full calculation with the one without FSI (dashed line). Since the s quark cannot be 
absorbed in medium the FSI can just increase the K + yield via processes like ttN — > K + Y 
(Y — E, A) for example. With decreasing formation time the primarily produced pions have 
a greater chance to produce K + in the FSI. As a consequence of this, a shorter formation 
time will lead to an increase of the nuclear K + photoproduction cross section as can be 
seen from the dotted line. An enhancement of the K + production cross section due to FSI 
cannot be explained by purely absorptive FSI as in simple Glauber theory where it would 
necessarily be interpreted as being due to a longer formation time of the K + . 



IV. SUMMARY 

High energy photoproduction off nuclei offers a promising possibility to study the physics 
of hadron formation. Necessary for such studies is a reliable model of the FSI to extract 
the formation time from the production cross sections. Whereas Glauber models allow for 
a straightforward implementation of the nuclear shadowing effect they usually have the 
disadvantage of a purely absorptive treatment of the FSI. As we have shown the latter can 
lead to a wrong estimate of the formation time. A more realistic treatment of the FSI is 
possible within a coupled channel transport model. We have presented a method to combine 
such an incoherent treatment of FSI with coherence length effects in the entrance channel 
which can easily be extended to higher energies and virtual photons. We have shown that in 
particular the production of mesons with long mean free path will be affected by the coupled 
channel effects in the FSI. 
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FIGURES 




FIG. 1. The two amplitudes of order a em that contribute to incoherent meson photoproduction 
in simple Glauber theory. The left amplitude alone would lead to an unshadowed cross section. 
Its interference with the right amplitude gives rise to shadowing. 
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FIG. 2. The nuclear photoabsorption cross section divided by Aa^N plotted versus the photon 
energy E~. The solid line represents the result of Reference [9] for the total photon nucleus cross 
section. The dashed line shows the contribution from incoherent reactions to a l °\ and is calculated 
using Eq. (2.12). More than 90% of the difference between <r^f and cr*^ is due to coherent p° 
photoproduction. 
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FIG. 3. The number density of nucleons that react with the <\> component (left side) and p° 
component (right side) of a 20 GeV photon for 208 Pb calculated using Eq. (2.14) . In both cases the 
nucleons on the front side of the nucleus shadow the downstream nucleons. This effect is stronger 
for the p° component because of its larger nucleonic cross section. 
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FIG. 4. Mass differential cross section for exclusive p production on 08 Pb at E y = 7 GeV. 
The meaning of the different curves is explained in detail in the text. All the curves, except the 
one with the explicitly given formation time ry = 0.8 fm/c have been calculated with r/ = 0. The 
fluctuations in the dashed and dash-dot-dotted curves are statistical only. 




|t| [GeV 2 ] 



FIG. 5. Calculated ^ for exclusive p° production on 208 Pb at E~ t = 7 GeV. The meaning of 
the different curves is the same as in Fig. 4. The structures in the curves at large \t\ are statistical 
only. 
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FIG. 6. Photoproduction cross section for K + on 208 Pb plotted as a function of the photon 
energy. The solid line represents the full calculation. The dash-dotted line shows the result without 
shadowing of the incoming photon, the dashed line the result without FSI and the dotted line the 
calculation without formation time. 
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